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SECTION  I 


INTRODUCTION 

A  study  recently  carried  out  by  Calspan  evaluated  the  characteristics 
of  a  number  of  radar  systems  with  the  possibility  that  such  a  system  might  be 
used  to  observe  satellites.  Two  important  characteristics  of  the  candidate 
systems  studied  for  this  application  were  that  they  were  narrowband  and  frequency- 
agile.  However,  no  data  were  available  on  satellites  that  could  be  used  in  simu¬ 
lations  of  these  systems.  Calspan  proposed  that  the  required  narrowband  data 
could  be  obtained  by  filtering  wideband  data  collected  by  the  Long  Range  Imaging 
Radar  (LRIR)  (also  known  as  the  Haystack  radar)  operated  by  Lincoln  Laboratory. 

In  this  report,  narrowband  means  bandwidth  such  that  individual  scatterers  on 
the  satellite  cannot  be  resolved. 

It  was  determined  that  the  work  that  needed  to  be  done  could  be  sep¬ 
arated  into  two  parts:  Phase  1  would  be  the  development  of  processing  software 
and  a  proof  of  the  concept  using  the  single  frequency  data  available  for  com¬ 
parison,  and  Phase  II  would  be  a  rather  detailed  analysis  of  frequency  agile 
signatures  corresponding  to  specific  candidate  radar  systems.  This  report  covers 
the  Phase  I  work  only,  and  due  to  a  change  in  requirements  for  this  data.  Phase  II 
will  probably  not  be  performed.  Section  2  of  this  report  is  a  mathematical  justi¬ 
fication  of  the  method,  with  some  details  given  in  Appendix  A.  Section  3  briefly 
describes  the  method  used  to  derive  the  required  signatures  from  the  data  avail¬ 
able  and  discusses  the  computer  software,  with  details  supplied  in  Appendix  B. 
Finally,  the  results  and  conclusions  are  given  in  Section  4,  and  may  be  summarized 
the  method  is  sound,  both  in  theory  and  practice. 


SECTION  II 


MATHEMATICAL  JUSTIFICATION  OF  APPROACH 

In  the  frequency  domain  a  wideband  radar  pulse  may  be  thought  of  as 
several  adjacent  narrowband  pulses.  Thus  it  is  quite  reasonable  to  expect  that 
individual  narrowband  pulses  could  be  obtained  from  wideband  pulses  by  some 
form  of  data  processing.  The  purpose  of  this  section  is  to  show  how  this  can 
be  done  and  to  define  the  restrictions  and  limitations  of  the  process,  both  in 
general  and  for  the  specific  case  of  the  LRIR. 


General  Approach 


Narrowband  radar  returns  can  be  derived  from  wideband  returns  by  a 
simple  filtering  operation.  This  can  be  shown  quite  easily  by  investigation 
of  the  received  signal  frequency  spectrum.  This  analysis  is  carried  out  in 
detail  in  Appendix  A  and  is  summarized  here. 


Consider  a  transmitted  waveform  with  transmitted  energy  E^  and 
spectrum'^  SQ(f)  reflected  from  a  point  target  as  in  Figure  2.1. a.  The 
received  signal  has  a  spectrum  B  SQ(f)  where  3  is  an  attenuation  factor. 

(For  simplicity,  Doppler  shift  has  been  omitted  from  this  discussion  although  it 
could  be  included  with  no  loss  of  generality.)  This  signal  together  with  additive 
white  Gaussian  noise  with  power  spectral  density  h.  Nq  is  mixed  down  to  the 
radar's  intermediate  frequency  and  input  to  a  filter  matched  to  the  Doppler 
shifted  signal.  The  filter  signal  output  spectrum  is  3 V  E^  [ SQ ( f ) | 2  with  an 
additive  noise  component  having  power  spectrum  ^  N^  | Sq (f ) ] 2 .  Note  that  here 
the  signal  (which  is  deterministic)  is  represented  by  its  spectrum  while  the 
noise  (which  is  a  stochastic  process)  is  represented  by  its  power  spectrum. 


The  output  signal-to-noise  ratio  (SNR)  is  given  by: 


SNR  = 


(Peak  Signal) 2 
Noise  Variance 


2  3  2E. 


Now  a  second  filter  H(f)  is  placed  after  the  matched  filter.  The  output  from 
the  second  filter  has  a  signal  component  spectrum  3  \  Ey  |Sg(f)|2  H(f)  and 
noise  power  spectrum  N^  | Sq ( f ) [ 2  | H ( f ) | 2 .  The  output  signal  components  for 
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Figure  1.  COMPARISON  OF  WIDEBAND  AND  NARROWBAND  RADAR  SYSTEMS 


this  system  (Figure  1-a)  are  »f¥o 


where 


SQ(f)  =  S0(f)  for  |f-f0|  <  B/2 

and  |  f+f ^ |  <_  B/2 
=  0  otherwise 

Next,  consider  the  system  shown  in  Figure  1-b.  The  wideband  signal 
with  spectrum  SQ(f)  is  first  passed  through  the  filter  H(f)  and  the  narrowband 
signal  "Y  E^,H(f)SQ(f)  is  transmitted.  The  received  signal  is  bV  E^H(f)S0(f) 
which  is  input  to  a  filter  matched  to  the  transmitted  signal.  The  output  has 
a  signal  component  with  a  spectrum  E^,|H(f)|2  |Sp(f)[2  and  a  noise  component 
with  power  spectrum  %  NQ  |.H(f)|2  |SQ(f)|2. 


The  noise  power  spectra  for  the  two  systems  shown  in  Figure  1  are 
identical,  but  their  signal  component  spectra  are  slightly  different.  The 
signal  components  will  only  agree  if  H(f)  =  | H (f ) [ 2  for  all  f.  Hence  if  11(f) 
is  defined  as: 


H(f)  =  1  for  |f-f0l  1  B/2 

and  |f+fgl  £  B/2 

=  0  otherwise 

Then  the  output  of  the  wideband  system  with  filtering  after  receiving  is 
identical  to  the  output  of  the  system  in  which  a  narrowband  signal  was  trans¬ 
mitted  and  received. 

Notice  that  the  output  signal-to-noise  ratio  for  either  system  is 

given  by: 

S  (f)[2df 

— - —  ] 

S0(f)|2df 

where  the  term  in  brackets  is  the  ratio  of  the  narrowband  transmitted  energy  to 
the  wideband  transmitted  energy. 


2  S2  Et  f, 

SNR  =  [  — 


N 


0 


l\ 
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The  conclusions  are  that: 


1.  Simulated  narrowband  data  can  be  extracted  from  wideband  field 
data  by  a  filtering  operation. 

2.  A  narrowband  transmitted  waveform  can  be  simulated  whose 
frequency  components  within  the  simulated  bandwidth  are 
identical  to  those  of  the  wideband  signal  within  the  same 
frequency  limits  and  zero  outside  these  limits. 

3.  There  is  a  loss  of  signal -to-noise  ratio  in  filtering  wideband 
data  to  obtain  narrowband  data  because  the  equivalent  transmitted 
energy  is  less  than  the  actual  wideband  transmitted  energy. 

2 . 2  Special  Considerations  for  Haystack  Radar 

In  the  general  discussion  given  above  it  is  assumed  that  the  wideband 
radar  return  passes  through  a  filter  exactly  matched  to  a  wideband  pulse.  The 
Haystack  radar  transmits  a  chirped  pulse  and  mixes  the  received  pulse  with  a 
replica  of  the  transmitted  pulse  delayed  in  time  to  correspond  to  a  return  from 
a  reference  range  and  offset  in  frequency  by  the  intermediate  frequency.  A 
return  from  the  reference  range  will  then  give  a  signal  at  the  center  of  the 
radar's  intermediate  frequency  band.  Returns  from  slightly  different  ranges 
will  result  in  constant  frequency  signals  with  frequency  proportional  to  range. 
However,  these  signals  are  also  slightly  displaced  in  time  from  the  reference 
signal.  The  receiver  filter  is  perfectly  matched  to  the  return  from  a  scatterer 
at  the  reference  range  at  a  particular  time,  but  at  that  time  it  is  not  perfectly 
matched  to  returns  from  scatterers  at  slightly  different  ranges.  The  purpose  of 
this  section  is  to  demonstrate  that  this  problem  is  negligibly  small  for  satel¬ 
lites  with  dimensions  of  interest  in  this  study. 
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Consider  the  way  in  which  the  chirp  pulses  are  collapsed  by  Haystack. 


Body  Center  at  Range  R^ 

Scatterer  i  at  range  R^  is  stationary  and  has  RCS 


Center  Frequency  =  f 
Swept  Bandwidth  =  B 
Pulse  Duration  =  T 


0 


We  use  complex  representation  to  simplify  equations.  The  transmitted 
signal  is  represented  as: 


s(t) 


j  27T[fQt  +  t2] 


for  -  T/2  <_  t  <  1/2 


.th 


The  return  from  the  i  scatterer  is 


'i(t)  =  VT~ 


B 


j2TT[f0(t-T.)  + 


CD 


(2) 


for  t ^  -  T/2  <_  t  +  1/2 


where  r.  =  2R-/c. 

l  i 


This  signal  is  beat  against  a  replica  of  s(t)  delayed  by  =  2Rq/c 
and  centered  at  a  lower  frequency  (f^  -  f jp) • 


=  ri^^  e 

■ 


-j2?r  [(fQ-f  jpKt-TQ)  +  2T^'t0^2^ 


(3) 


" j 27Tf ip  T0  j27r[f0(x0-Ti)  +  (Ti2-T02^  J 2w tf ip  *  T  (Ti_T0)]  1 

e  e  ^ 


for  -  T/2  +  t.  <  t  <  T/2  +  r. 

l  —  —  l 
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Hence  a  scatterer  at  range  i  yields  a  CW  pulse  centered  at  t.  at  frequency 

B  ^ 

fjF  -  j  (ty  -  x0).  Note  that  the  duration  of  the  chirp  affects  only  the  limits 

of  the  return;  the  phase  and  frequency  of  the  return  are  fully  determined  by  the 

chirp  rate  (^) • 

The  return  from  several  scatterers  comprising  a  realistic  target  can 
be  found  by  summing  the  y^(t)  for  each  of  the  scatterers. 


Figure  2  depicts  the  individual  components,  y^Ct),  (and  their  sum¬ 
mation,  y(t))  of  the  return  from  a  scatterer  located  at  t^,  the  reference  range, 
one  located  slightly  closer  to  the  radar  at  t  and  one  located  slightly  farther 
from  the  radar  at  t+ .  The  figure  shows  clearly  the  effect  of  the  slight  dis¬ 
placement  of  the  individual  scatterers  from  each  other.  Let  us  define  a  segment 
of  y(t)  centered  at  tQ  with  duration  6  : 


y'(t)  =  y(t) 

=  0  otherwise 


f°r  T0  -  I  i  tQ  +  T0  -  \  <  t  <  tQ  ♦  T0  +  |  <  T0  +  J  (4) 


.th 


y!(t)  =  Va.  e 


The  contribution  of  the  i  scatterer  (i=+,  0,  or  -)  in  this  signal  is: 
j2TTfIF  T0  j27T[f0(T0-Ti)  +  ^  (T^-To2)]  j  2  TT  [  f  j  p  -  £  t 


(5) 


where  tQ  +  tq  -  ^  <  t  <  t„  +  t0  +  5y 


-j2^fIF  tq  j2it[f0(Xo-Ti)  +  ^  52nfiF-  f  (^i-t0)][t+t „] 

y'  (t+tQ)  = V  oi  e  e  e  (6) 


■v^r 


”-’2lTflF^T0't(P  j2lT[(fo+T  t0^T0_Ti-*  +  7f  ^Ti2_T02^  j2lTtfIF'  T^Ti"T0^t 
ai  e  e  e 

(7) 


where  Tq 


±T0 


+ 


5_ 

2 
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Figure  2  COMPONENTS  AND  COMPOSITE  OF  CHIRPED  RETURN  FROM  3-SCATTERER  TARGET 


Equation  7  has  the  same  form  as  the  completely  general  Equation  3. 

This  shows  that  except  for  fringe  or  end  effects,  the  signal  from  a  multiple 
scatterer  target  is  in  all  other  respects  completely  general.  Therefore,  when 
the  chirped  radar  pulse  is  very  long  in  comparison  to  the  range  extent  of  the 
target  (At<<5),  as  it  is  for  the  LRIR  looking  at  satellites  of  reasonable  size, 
the  filtering  used  with  wideband  radar  returns  is  negligibly  different  from  a 
true  matched  filter  and,  therefore,  should  not  have  a  significant  effect  on  the 
derivation  of  narrowband  signals  from  the  wideband  data. 

SECTION  III 
DESCRIPTION  OF  METHOD 

As  the  discussion  in  the  previous  section  implies,  filtering  wideband 
radar  returns  to  obtain  narrowband  returns  is  a  relatively  simple  process.  Most 
of  the  effort  in  this  study  was  related  to  the  form  in  which  the  data  were  avail¬ 
able.  In  order  to  minimize  the  effort  required  to  obtain  the  data,  understand 
and  compensate  for  radar  dependent  characteristics,  remove  translational  motion 
of  the  satellite,  and  be  certain  of  the  validity  of  the  data,  the  most  practical 
form  for  transferring  the  wideband  data  from  Lincoln  Laboratory  to  Calspan  was 
in  the  form  of  two  dimensional  images.  Since  these  images  are  the  result  of 
reversible  processing  of  the  needed  wideband  data,  much  of  the  software  developed 
in  this  study  was  aimed  at  reversing  many  of  the  processing  steps  involved  in 
creating  the  images. 

One  way  to  describe  the  methods  used  for  this  processing  is  briefly  to 
review  this  imaging  process.  Figure  3  shows  the  ramp  waveform  transmitted  by 
by  LRIR.  The  received  waveform  is  mixed  with  a  delayed  (to  an  appropriate  ref¬ 
erence  range)  replica  of  the  transmitted  waveform  to  generate  the  uncollapsed 
pulse.  This  pulse  is  digitized,  Doppler  corrected,  weighted  and  Fourier  trans¬ 
formed  to  obtain  a  high  resolution  pulse,  where  successive  samples  from  the 
transform  represent  range  samples  near  the  reference  range.  Samples  of  many 
pulses  at  the  same  range  relative  to  a  reference  point  fixed  on  the  target  are 
weighted  and  Fourier  transformed  in  the  cross-range  dimension  to  obtain  Doppler 
or  cross-range  resolution  of  target  scatterers.  This  final  step  is  repeated 
for  each  sampled  range,  and  the  composite  represents  a  two-dimensional  image  in 
the  range-cross-range  plane. 
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CROSS-RANGE  AMPLITUTE  AMPLITUDE  TIME 


RAMP  WAVEFORM 


rmw 


UNCOLLAPSED  PULSE 


FREQUENCY 


HIGH  RESOLUTION  PULSE 


RANGE 


RADAR  IMAGE 


Figure  3  REVIEW  OF  IMAGING  PROCESS 
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The  software  implemented  in  this  study  reversed  most  of  the  steps  in 
image  processing.  It  is  important  to  note  that  two  rather  involved  steps  did 
not  have  to  be  reversed  as  they  would  have  been  required  even  if  the  raw  data 
had  been  directly  available  to  Calspan.  The  first  of  these  was  the  compensation 
discussed  in  Appendix  A  to  remove  the  apparent  range  shift  due  to  Doppler  fre¬ 
quency.  The  second  was  the  very  accurate  range  tracking  of  the  target  that  was 
required  to  ensure  that  the  same  frequency  could  be  selected  from  successive 
radar  pulses. 

The  software  implemented  is  summarized  by  the  block  diagram  in  Figure  4. 
The  plotting  of  the  two-dimensional  images  was  done  to  ensure  that  the  data  were 
being  read  properly  from  the  magnetic  tapes.  The  narrowband  images  were  plotted 
to  facilitate  eyeball  correlations  among  signatures.  The  majority  of  the  soft¬ 
ware  performed  the  reverse  image  processing  steps  described  above  and  presented 
various  outputs  to  permit  checking  of  intermediate  results.  Simple  statistics 
were  calculated  to  permit  comparisons  among  actual  narrowband  signatures  and 
derived  narrowband  signatures.  This  software  is  described  in  detail  in  Appendix  B. 

SECTION  IV 

RESULTS  AND  CONCLUSIONS 

The  primary  purpose  of  this  investigation  was  to  demonstrate  that 
narrowband  signatures  derived  from  wideband  data  are  a  good  representation  of 
actual  narrowband  signatures  in  practice  as  well  as  in  theory.  This  section 
describes  the  results  of  performing  the  required  comparisons. 

The  nature  of  the  LRIR  precluded,  at  reasonable  cost,  obtaining  wide¬ 
band  and  actual  narrowband  data  simultaneously.  The  data  available  for  this 
investigation  are  summarized  in  Table  1.  Since  the  geometry  differences  from 
pass  to  pass  would  make  the  signatures  quite  different,  only  a  qualitative  pattern 
lobe  correlation  could  be  made  between  different  passes;  quantitative  comparisons 
could  be  made  on  a  statistical  basis.  Therefore,  it  was  decided  to  concentrate 
on  satellite  10352  since  there  were  three  actual  narrowband  signatures  available 
to  compare  to  each  other  and  to  the  signature  at  the  same  frequency  that  could 
be  derived  from  wideband  data. 
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NARROWBAND 


SIGNATURE 
FROM  HAYSTACK 


Figure  4  DATA  REDUCTION  SOFTWARE  FUNCTIONAL  DIAGRAM 


TABLE 


TOTAL  124  IMAGES 


The  data  from  the  three  narrowband  passes  had  been  collected  at  either 
100  or  200  pulses  per  second  and  a  frequency  of  10.0  GHz.  Since  this  was  much 
higher  than  the  effective  pulse  rate  (three  to  four  pulses  per  second)  used 
in  the  imaging  process,  groups  of  narrowband  data  were  averaged  to  obtain  a 
comparably  low  sample  rate.  This  averaging  also  eliminated  a  non-random  noise 
that,  Lincoln  Laboratory  personnel  agreed,  was  a  figment  of  the  radar  operating 
in  the  narrowband  mode  rather  than  a  characteristic  of  the  satellite.  These 
data  were  then  plotted  and  statistics  calculated  to  permit  comparison  to  each 
other  and  to  the  signatures  derived  from  the  wideband  data. 

The  wideband  data  from  the  pass  of  Object  10352  were  processed  as 
described  in  Section  3,  selecting  the  same  10.0  GHz  center  frequency  as  used 
above.  The  signature  segments  derived  from  each  image  were  concatenated  to 
obtain  a  signature  comparable  to  that  from  a  narrowband  pass.  This  signature 
was  plotted  and  its  statistics  calculated. 

The  plots  from  the  three  actual  narrowband  passes  and  the  narrowband 
data  derived  from  the  wideband  pass  were  compared  by  "eyeball"  correlation  to 
see  if  they  were  generally  similar.  It  was  found  that  all  the  signatures  had 
generally  similar  characteristics.  Figures  5  and  6  show  segments  of  two 
of  these  signatures  The  similarities  and  differences  are  quite  reasonable 
for  two  different  passes  of  the  same  satellite  as  seen  at  a  radar  frequency 
of  10.0  GHz.  One  other  characteristic  that  could  be  seen  but  not  quantified 
was  the  match-up  of  the  segments  of  the  derived  signature  at  the  ends  of  the 
individual  images.  The  18  images  gave  17  junctions  of  which  three  had  a  rela¬ 
tively  long  gap  between  the  end  of  one  image's  data  and  the  beginning  of  the 
next,  eleven  had  short  gaps  and  three  had  overlapping  data.  In  the  cases  of 
short  gaps  or  overlaps,  the  excellent  continuity  and  retracing  of  the  sig¬ 
nature  plot  was  an  indication  that  the  signature  segments  were  independent 
of  the  individual  images  and  that  there  were  no  significant  effects  at  the 
ends  of  the  segments  of  the  signature.  Although  these  two  observations  did 
not  directly  validate  the  derived  signature,  they  did  indirectly  add  credence 
in  a  double  negative  sense:  the  signature  did  not  show  a  problem. 
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Figure  5  PORTION  OF  NARROWBAND  SIGNATURE  OF  OBJECT  10352  AT  10.0  GHz 


Analysis  of  the  signature  statistics  led  to  a  more  positive  indication 
of  the  validity  of  the  derived  signature.  Table  2  shows  several  statistics 
that  compare  the  derived  signature  to  the  three  actual  signatures.  Clearly,  the 
derived  signature  fits  well  with  the  actual  signatures,  considering  the  varia¬ 
tion  among  the  actual  signatures.  The  other  statistic  available  for  comparison 
was  the  distribution  of  RCS  values  over  the  signatures.  The  normalized  dis¬ 
tributions  are  shown  in  Figure  7,  and  the  corresponding  cumulative  distributions 
are  shown  in  Figure  8.  Again,  the  derived  signature  fits  well  within  the 
bounds  established  by  comparison  of  the  actual  narrowband  signatures. 


TABLE  2 


COMPARISON 

OF  ACTUAL 

AND  DERIVED  RCS 

SIGNATURES 

Pass 

Mean 

RCS 

Standard 
Deviation 
of  RCS 

Max imum 

RCS 

Minimum 

RCS 

Actual  Narrowband 

11/8/78 

1.3 

5.7 

12 

-25 

Actual  Narrowband 

11/9/78 

-1.5 

5.9 

11 

-27 

Actual  Narrowband 

11/8/78 

0.9 

5 . 6 

12 

-23 

Derived  from  Wideband 
11/7/78 

1.4 

5.4 

13 

-25 

The  major  conclusion  that  can  be  drawn  from  this  investigation  is 
that,  theoretically  and  practically,  narrowband  signatures  derived  from  wide¬ 
band  data  are  a  very  good  representation  of  actual  narrowband  signatures.  The 
most  important  restriction  is  that  the  signal -to-noise  ratio  (SNR)  of  the  de¬ 
rived  signature  is  lower  than  the  original  wideband  SNR  by  the  ratio  of  the 
bandwidths,  if  the  wideband  signal  spectrum  is  flat,  as  is  typical.  There  is 
no  restriction  on  the  complexity  or  number  of  radar  scatterers  of  the  target. 
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RCS  DISTRIBUTION 
10  GHz  SIGNATURE,  OBJECT  10352 


Figure  7  RCS  DISTRIBUTION 


CUMULATIVE  RCS  DISTRIBUTION 
10  GHz  SIGNATURE,  OBJECT  10352 


Figure  8  CUMULATIVE  RCS  DISTRIBUTION 


The  value  of  this  capability  to  derive  narrowband  information  from 
wideband  data  is  multifold.  1)  Narrowband  data  can  be  obtained  without  the 
necessity  of  modifying  a  wideband  radar  transmitter  and  receiver.  2)  Narrow- 
band  data  at  multiple  frequencies  can  be  obtained  simultaneously  from  the  wide¬ 
band  data,  permitting  different  radars  and/or  frequency-agile  radars  to  be  com¬ 
pared  under  identical  geometric  conditions  and  with  zeal  rather  than  simulated 
data.  An  example  of  a  return  at  a  second  frequency  is  given  in  Figure  9; 
the  sample  corresponds  to  that  shown  in  Figure  6.  3)  Although  not  specifically 

discussed  here,  phase  coherence  is  not  lost  in  this  filtering  process,  which 
means  that  coherent  waveforms  can  be  evaluated  using  data  derived  by  this  process. 


PORTION  OF  DERIVED  NARROWBAND  SIGNATURE  OF  OBJECT  10352  AT  10.032  GHz 


APPENDIX  A 


DERIVATION  OF  NARROWBAND  SIGNALS  FROM  WIDEBAND  SIGNALS 


This  appendix  shows  that  narrowband  radar  returns  can  be  derived  from 
wideband  returns  by  a  simple  filtering  operation.  This  is  done  by  showing  that 
the  output  from  filtering  a  wideband  return  is  identical  to  the  output  that 
would  have  been  obtained  if  the  radar  had  transmitted  and  received  a  narrowband 
signal. 


A  wideband  signal  S(t)  is  transmitted;  it  has  the  form: 
S (t)  =  V  2  ET  g(t)  cos  [2nf0t  +  4>(t)  ] 

where  /“  f2(t)  dt  =  1 


g(t) 


and 


transmitted  signal  energy 
amplitude  modulation  function 
phase  modulation  function 
RF  frequency 


The  signal  is  reflected  from  a  point  target  and  the  return  can  be 
represented  as: 

r (t)  =  B  VTet  fc(t)  cos  [2Mf0+  f„)  t  +  ♦(t)l  +  n(t) 

where  6  =  channel  attenuation 


fp  =  Doppler  frequency  shift 

n(t)  =  zero  mean  white  Gaussian  noise  with  spectral  density  h  NQ 


Here  it  is  assumed  that  the  target  moves  no  more  than  a  small  fraction 
of  a  range  resolution  cell  during  the  coherent  duration  of  the  waveform.  If  this 


1 

H 


is  not  true,  as  with  the  Haystack  radar  looking  at  satellites,  a  simple  correc¬ 
tion  is  made  to  compensate  for  the  known  Doppler  shift. 


The  frequency  spectrum  (amplitude  only)  of  the  signal  portion  of 
this  return,  S(f),  is  shown  in  Figure  A1  (positive  frequency  only)  with  the  noise 
power  spectral  density,  Pg^ff) »  shown  below  it.  The  combined  signal  plus  noise 
return  is  supplied  to  a  matched  filter  in  the  receiver  which  effectively  multi¬ 
plies  the  frequency  spectrum  of  the  signal  SQ(f)  by  S*(f)  and  the  noise  power 
spectrum  density  is  multiplied  by  S^ff) 


where : 


Vf+V 


/“g(t)  cos 


-j2Tlft 

[  2ti  (f  0+fD)  t  +  <)>(t)]  e  dt 


The  spectrum  at  the  output  of  the  matched  filter  may  be  considered  as 
the  combination  of  signal  and  noise  related  terms: 

|Sg(f)|z  and  respectively,  and  are  shown  in  Figure  A2. 


Now  assume  that  this  signal  is  next  applied  to  an  ideal  filter  with 
center  frequency  f^  and  bandwidth  11.  The  signal  spectrum  portion  of  the  output 
of  this  filter  is: 

|SQ(f)|2  =  | S(f ) | 2  for  |f+fc|  <  B/2 

and  |  f-f  |  _<  B/2 
=  0  otherwise 

and  the  noise  power  spectrum  portion  is: 

Wf>  ’  >Wf>  •  X  lS<f>l2  f°r  lf*fol  iB/2 

and  |  f-f  |  <_  B/2 


0  otherwise 


This  return  from  the  ideal  filter  is  equivalent  to  the  output  of  the 
matched  filter  if  the  input  signal  was  originally 

f0+B/2  j2irft 

S(t)  =  2  Re  {  /  Sft(f)  e  df  } 

f0-B/2  0 

It  is  exactly  the  return  that  would  have  resulted  if  a  narrowband  signal  (fre¬ 
quencies  from  f^-B/2  to  fQ+B/2)  had  been  transmitted  and  received.  In  other 
words,  filtering  a  wideband  return  after  receiving  will  obtain  the  same  result 
as  transmitting  and  receiving  a  narrowband  signal. 
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APPENDIX  B 


SOFTWARE  DESCRIPTION  AND  USERS' GUIDE 


I .  Physical  Characteristics 


Tape  #AXDSRC  is  a  9-track,  non-labeled,  800  BPI  tape  with  the  following 
DCB  information:  RECFM=FB,  LRECL=80,  and  BLKSIZE=8000.  It  is  IBM  compatible  and 
contains  3  files. 


II.  Format 


Each  record  on  the  tape  is  a  card  image  of  80  characters  and  may  be 
read  in  the  same  manner  used  to  read  a  deck  of  cards. 


III.  Contents 


The  content  of  each  file  is  as  follows: 

1.  File  1 

This  file  contains  the  software  used  to  derive  the  narrowband  data 
from  the  wideband  coherent  images.  It  includes  a  main  driver,  a  block  data  routine, 
processing  routines  and  output  routines.  The  prologues  (including  inputs  and 
outputs)  are  given  in  Table  B1 .  Below  is  a  brief  description  of  each: 

A)  MATH  -  reads  inputs  and  controls  program  flow  through  the 

processing  and  output  routines. 

B)  BLOCK  -  a  block  data  to  initialize  constants  in  common  /LORINF/ 

for  ladder  plots. 

C)  Processing  routines: 

INV1  -  performs  an  inverse  Fourier  transform  on  complex 

image  data  in  the  cross-range  direction  to  generate 
wideband  pulse  shapes. 
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INV2 


performs  an  inverse  Fourier  transform  on  the  complex 
pulse  shapes  to  generate  frequency  spectra. 

KAISER  -  generates  a  set  of  Kaiser  weights  for  the  frequency 
spectra. 

PRDIST  -  calculates  the  mean,  standard  deviation,  minimum, 
maximum,  frequencies  of  occurrence  and  cumulative 
distribution  of  a  set  of  data. 

D)  Output  Routines: 

LORINI,  LDRPLT,  PLT  -  initializes  for  and  generates  the  side¬ 
band  ladder  plots. 

MAGPLT  -  generates  a  printer  plot  of  the  magnitudes  of  a 
2-dimensional  image. 

0UT1  -  generates  a  listing  of  data  statistics  calculated  in 
subroutine  PRDIST. 

E)  The  following  are  routines  referenced  in  this  program  but  not 
included  on  this  tape  due  to  particular  machine  dependency: 

FOURT  -  performs  Cooley-Tukey  Fast  Fourier  Transform 

CORE  -  allows  use  of  I/O  statements  to  perform  conversions 
on  a  variable  list  into  core  instead  of  to/from  an 
I/O  device 

CTIME2,  DATE,  JBNAME  -  returns  time-of-day,  date  and  name  of 

job  into  core 

PLOT,  SYMBOL,  NUMBER,  EFPLOT  -  Calspan' basic  plot  software 

PLTTME  -  see  File  3. 

2.  File  2 

This  file  contains  the  software  used  to  convert  raw  wideband  radar 
data  into  2-dimensional  contour  images.  The  software  is  composed  of  a  main 


driver/processing  routine  and  a  contour  plot  subroutine.  The  prologues  de¬ 
scribing  each  routine  and  its  usage  are  given  in  Table  B2.  Below  is  a  brief 
description  of  each: 

A)  IMAG2  -  reads  input  data  from  tape  and  cards,  processes  the 

data  for  plotting,  calls  contour  subroutine  and  then 
writes  identifying  information  on  contour  image. 

B)  MCCNTR  -  the  contour  plot  subroutine,  employs  Calspan's  basic 

plot  software  (PLOT,  AXIS,  SYMBOL,  NUMBER)  to  generate 
the  2-dimensional  contour  images. 


3.  File  3 

This  file  contains  the  software  used  to  process  narrowband  radar  data 
and  output  results  in  the  form  of  plots  and  statistical  summaries.  The  software 
is  composed  of  a  main  driver,  processing  routines  and  output  routines.  A  brief 
description  of  each  routine  will  be  given  below,  while  the  prologues  describing 
each  program  can  be  found  in  Table  B3. 

A)  NBAND  -  reads  parameter  cards,  advances  tape  to  initial  position, 

controls  flow  through  processing  and  output  routines. 

B)  Processing  Routines: 

SMOOTH  -  reads  narrowband  tape  of  RCS  data  in  dBsm  units, 

converts  to  square  meters  and  calculates  average  RCS 
for  the  user  selected  number  of  data  points. 

CONVRT  -  converts  an  array  of  RCS  data  points  from  square 
meters  to  dBsm. 

PRDIST  -  calculate  the  mean,  standard  deviation,  minimum, 
maximum,  and  distribution  of  a  set  of  data. 

C)  Output  Routines: 

0UT1  -  outputs  to  the  printer  the  statistics  calculated  in  PRDIST 

PLTTIME  -  outputs  narrowband  plot  of  RCS  vs  time. 
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PLTASP  -  outputs  narrowband  plot  of  RCS  vs  aspect  angle  or 
projected  aspect  angle. 

Note:  As  pointed  out  in  File  1  description,  the  programs  on  both 
File  2  and  File  3  call  for  Calspan's  basic  plot  package  composed  of  PLOT,  AXIS, 
SYMBOL,  NUMBER, --which  have  not  been  included  on  this  tape  due  to  particular 
machine  dependency. 
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TABLE  B1 


PROLOGUES  TO  ROUTINES  ON  FILE  1 


cc** 

cc* 

cc * 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 


*****  **  ****  ****************  MAIN  ******************♦**************< 


PURPHSfc 

ro  DERIVE  N.B  DA  T  A  FROM  WH  DATA  THROUGH  INVERSE 
h'UHM'j  . 


:OUKlER  TWANS-*i 


M  E  THUD  *» 

AN  INVERSE  r'UUWlEK  FKANSr'UKM  IS  PERFuRMED  UN  THl  COMPLEX  *» 

image  in  the  CROSS-RANGE  DIRECTION  FOR  EACH  range  GATE  *1 

WeOUeSTeU.  THIS  GENERATES  T  ME  ORIGINAL  WB  MJLSt  SHAPES.  AN  *< 

INVERSE  FFT  IS  THtN  OONl  UN  lACH  PULSE  SHAPE  T'l  GLNeRATL  *» 

FREQUENCY  SPECTRA.  VARIOUS  TAPE*  LIST  AND  PLOTS  ARE  PRODUCED  *i 
AT  EACH  STAGE  OF  THE  PROCESSING.  FINALLY*  STATISTICS  ARE  DO^*< 
ON  A  SINGLE  FNLOUENCY  CIFSTAT).  *t 

*i 

PROGRAM  INPUTS:  *< 

tapes  -  UNIT  l:  WH  COMPLEX  IMAGE  TAPE  FROM  Ll  *‘ 

CARDS  -  UP  TO  3  NAMELISTS  ARE  READ  t  *< 

C  INPUT  l  -  REOUIRLO  *« 


ngatci 

ngatls 


FIRST  RANGE  GATE  IN  IMAGE  TO  BE  PROCESSED  *‘ 
-  NUMBER  OF  GATES  TO  BE  PROCESSED  (MUST  BE  *' 
DIVISIBLE  BY  4)  *( 

NIMAGS  -  NUMHER  OF  IMAGES  TO  PROCESS  *< 

SLL  -  SI  OELObt  LEVLL(DB)  FOR  WEIGHTING  *< 

MGPLTS—  LOG.  INDICATOR  FUR  MAGNITUDE  PLOTS  UK  T  rfc  *> 
IMAGE*  wB  PULSES  &  Fr£u(FJR  EACH  IMAGE)  *< 

4BLUK  -  LOGICAL  INDICATOR  FOR  WB  LADDER  PLOT  * 

wbtape-  lug.  indicator  for  wu  pulse  shape  tape.  *• 

WlJLST  -  LOG.  INDICATOR  FOR  D I  AG .  LIST  OF  1  PULSE  *> 

IPLlST-  PULSE  NUMBER  TO  BE  LISTED  *' 

ndplt  -  log.  indicator  for  nb  rcs  vs.  time  plot  *• 

NuTAPfc-  log.  indicator  for  no  frequency  tape  *• 

nblst  -  lug.  Indicator  fur  diag.  list  of  j  fred.  ** 
I FL  I  ST  —  FkED  NUMBLR  TO  HE  LlSTLD  *> 

IFPLOr-  FRED  NUMBER  TO  BE  PLOTTED  UN  RCS  PLOT  * 

IFSTAT-  FREU  NUMBLR  FOR  STATISTICS  *’ 

RCSL  -  LOW  BIN  RCS  VALUE  FOR  HISTOGRAM  *> 

RC  SR  -  HIGH  BIN  RCS  FUR  HIST.  *' 

rCSInc-  rcs  Increment  for  hist.  *• 

CALIB  -  CALIBRATION  FACTOR  FOR  FRED  DATA  *• 

ELDRIN  -  RLA.D  ONLY  IF  W3LDK=T  * 

XLEF  T  *  X  R I GHT  *  X  LEN  *  YBOT  »  Y  T  OP ♦ YLEN • LABRL  P • RUNGS  P  * 

**fur  description  of  these  see  common  cdrinf  beluw  *• 

ETMEPLT  -  READ  (BY  SUBROUTINE  PLTTMd)  ONLY  IF  NBPLT  =T*« 

RCSLNG-  LENGTH! IN)  OF  RCS(Y)  AXIS  * 

RCSSCL-  INCRLMi.NT  PER  INCH  (DO)  FOR  RCS  * 

TMFSCL-  I-NCREMtNT  PER  INCH(SCC)  FUR  TIMe(X)  AXIS  * 

T  APE  ID—  TITLE  TO  HE  PLACED  UN  RCS  PLOT  * 

RCSTTL-  TITLE  FOR  RCS  AXIS  * 

RCSMIN-  MIN.  VALUE  FOR  RC  S  * 

* 

PROGRAM  OUTPUTS:  * 

TAPlS  -  UNIT  ?:  wo  PULSE  SHAPES  (IF  WBTAOE-T)  * 

UNIT  3:  DERIVED  NB  DATA  FOR  ACL  FREQUENCIES  (IF  NOT APl' 

=  T  )  *i 

,?  PUSSIHLE  (NOT  BOTH  IN  1  RUN):  *< 

1)  WH  PULSE  SHAPE  LADDER  PLOT  (IF  WHLDR—  T  )  *' 

d)  NB  RCS  VS.  TIME  PLOT  FUR  FRED  *  IFPLOT  (IF  NbPLT  =T )  * 

LlST.>  -  PRINTER  PLUTS  OF  MAGNITUDE!  IF  MGPL  T  b=  T  >  OF:  *' 

VB  IMAGES*  Wt»  PULSES  AND  NB  FREQUENCIES  *' 

DIAG.  LIST  UF  PUlSl  *  IPL I  ST  (IF  WrfCST  =  T)  * 

DIAG.  CIST  OF  FRLU  *  IFLIST  (IF  N3lST  =  T)  * 

STATISTICAL  RESULTS  FOR  EACH  IMAGE  * 

* 

COMMuN  VARIABLES  USED  ♦ 

/LOHINF/  -  VARIABLES  USED  FOR  WB  LADDER  PLOT  * 

• 

REMARKS  * 

-  TJ  PRODUCE  A  Nil  LADDER  PLOT*  PL  UT  TE  R  =  CAL  C  :JMP  *  Ll)NG=V  ,  * 

PA^E«=SG  SHOULD  BE  SPECIFIED  ON  THE  EXEC  CARD  * 

-  TO  PRODUCE  A  NU  RCS  PLOT*  PLOTTER=C ALCUMP *LUNG  =  X *PAPER=SG  * 

SHOULD  BE  USED  * 

-  AT  THIS  TIME*  THE  PROGRAM  CAN  DO  EITHER  THE  LADDER  PLOT  OR  * 

T  Hi.  RCS  PLOT*  NOT  BOTH  * 


plots 


>  ^  W  ti 


>  o 


vT3' 

'id. 


,«*V 
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TABLE  B1 


PROLOGUES  TO  ROUTINES  ON  FILE  1  (CONT.) 


CC* 

CC* 

SUBROUT 

cc* 

CTI 

cc* 

K  A  I 

cc* 

INV 

cc* 

cc* 

cc* 

A  UtHOR/ 

cc* 

C.M 

cc* 

MAW 

cc* 

DISKL1B  •  lCMD,  f) I SKL IB* 


'.9  79 

cc**  ***  ************************  block  *******  **********************  **> 

cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 

CC**  **********************  *******  ************************************  ** 

cc***** ************************  INVI  ******************************  ♦* 

cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
ec* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 


block 
SUBROUTINE  BLOCK 

PURPOSE 

TO  INITIALIZE  VARIABLES  IN  COMMON 
COMMON  VARIABLES  OSLO 

VARIABLES  ARE  INITIALIZED  IN  THE  FOLLOWING  COMMONS: 

LOR  INF 

AUTHOR/PROGRAMMER 
C .M.SZCZESNV 
MAR  1979 


SUBROUTINE  INVI 

PURPOSE 

ro  PERFORM  AN  INVERSE  FOURIER  TRANSFORM  ON  COMPLEX  IMAGE 
OATA  IN  THE  CROSS-RANGE  DIRECTION,  TO  GENERATE  WB  PULSE 
SHAPES* 

CALLING  SEQUENCE 

CALL  INVI  IV.NF.NG.NTCJL  ) 

DESCRIPTION  OF  PARAME  TERS 

inputs: 

NF  -  NO*  OF  CRUSS-RANGE  CELLS  IN  THE  IMAGE 

NG  “  NO.  OF  RANGE  GATES  IN  THE  IMAGE 

N  TO L  “  ACTUAL  NO.  UF  PULSES  USED  TO  GENERATE  THE  IMAGE 

hybrids: 

V  -  COMPLEX  ARRAY  CONTAINING  THE  2-0 1  MENS IONAL  IMAGE 

COMMON  VARIABLES  USED 
NONE 

SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
FOURT  -  CALSUBS  FOURIER  TRANSFORM  ROUTINE 

A UTHUR/ PROGRAMMER 
M  .BURNS 
MAR  1979 


CC***** ******  ******************************************  ******?**■*****  ** 
CL«***« ************************  INV2  ******************************** 


CC* 

CC*  SUBROUTINE  I  NV2 

CC* 

CC*  PURPJSE 

cc*  TO  PERFORM  AN  INVERSE  FOURIER  TRANSFORM  ON  COMPLEX  IMAGE 

CC*  DATA  IN  THE  RANGE  DIRECTION  TO  GENERATE  FREQUENCY  SPECTRA 

CC*  FROM  THE  W B  PULSE  SHAPES. 

CC* 

cc*  CALLING  SEQUENCE 

CC*  CALL  INV2  <V.NF.NG#NTOL.*T2*CALlD J 

CC* 

CC*  DESCRIPTION  OF  PARAMETERS 

cc*  inputs: 

cc*  NF  -  NO.  OF  CROSS-RANGE  CELLS  IN  THE  I 

CC*  NG  -  NO.  OF  RANGE  GATES  IN  THE  IMAGE 

CC*  NTOL  “  ACTUAL  NO.  'IF  PULSES  USED  TO  GENE 

CC*  Wf2  -  ARRAY  OF  WEIGHTS  FOR  THE  RANGE  OI 

CC*  C  AL  18—CALlORAT  (  ON  FACTOR  FOR  SPECTRAL  DA 

cc*  hybrids: 

CC*  V  -  COMPLEX  ARRAY  CONTAINING  THE  2-01 

CC* 

CC*  COMMON  VARIABLES  USED 
CC*  NONE 

CC* 

CC*  SUBROUT  INE  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
CC*  FtlURT  -  CALSUBS  FOURIER  TRANSFORM  ROUTINE 


MAGE 

RATE  THE  IMAGL 
MENS  ION 

r  a 

MENS IONAL  IMAGE 
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cc*  * 

CC*  AUTHOR/ PROGRAMMER  * 

CC*  M  •BURNS  « 

CC*  MAR  1979  _  . .  .  .  _ 

cc*  “  *c 

CC*****  »»»»»» <t» »»»»»»»»»♦»»»»»»♦♦  [**»**<L*»*_**  ♦♦»»»♦  **»*»♦  +** 

cc *****  ********  *******  *********  KAISER  **********  ******  *********  ****** 

CC* 

CC* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 
cc* 

CC***** %***************************************$*********************  *** 
CC*****  ******** ***********4****  MHO I  ST  ******************************* 

cc*  * 

CC*  SUBROUTINE  PROIST  * 

CC*  * 

CC*  PURPOSE  * 

c**  ^VALUATION  of  data  statistics  * 

C**  1 >  MEAN  VALUE  * 

C**  <t)  STANDARD  DEVIATION  * 

C**  3)  MAXIMUM  * 

C**  4)  MINIMUM  * 

C**  5)  HISTOCRAMIFREUUENCIES  OF  OCCURRENCE)  * 

C**  ft)  CUMULATIVE  DISTRIBUTION  * 

CC*  * 

CC*  CALLING  SEQUENCE  * 

CC*  CALL  PRD  1ST  ( X  *NP *A*R*C*D*EX )  „  * 

CC* 

CC*  DESCRIPTION  of  parameters  * 

CC*  inputs:  • 

C**  X  IS  THL  NAME  FOR  THE  OATA  ARRAY  * 

cc*  np  is  the  number  of  data  points  * 

C**  A.H.C.....  DEFINE  CELL  SIZES  AS  * 

C**  A(B)C(D>E  -  <  THAT  IS*  A  TO  C  IN  INCREMENTS  OF  SIZE  B «  ♦ 

C**  THEN  TO  E  IN  INCREMENTS  OK  SUE  D)>  * 

C**  FOR  A  {  U  >C  *  SET  E=D=?C  OR  SET  Es  0*0*  • 

CC*  * 

CC*  COMMON  VARIABLES  USED  • 

cc*  outputs:  _  - 

CC*  /LCPMOB/  SCL  OESC.  OF  PARAMETERS  IN  REMARKS  * 

CC*  * 

CC*  REMARKS  * 

c**  results  are  left  in  labeled  common  lcp.rob.  defined  by  * 

c**  COMMON/LCPROB/  NBUCM  10  I  )  *D  I  ST (  1  01  >  #CELL1  (  1  02).  »XMAX  *XMI  N.  * 

C**  XBA«*DEV .Kl *NX  • 

C**  NBUCKI  1 )  IS  THE  NO.  OF  OCCURRENCES < I  NT  EGER )  FOR  THE  ITH  CELL*  • 

C**  DISTU)  IS  THE  CUMULA  TIVE  DI  S TRIBUTE ON (FLOATING  POINT  VALUE  IN* 

C**  INTERVAL  0.0  TO  1.0)  CORRESPONDING  TO  THE  RIGHT  HAND  * 

C**  EDGE  OF  THE  ITH  CELL.  * 

C**  CELL  HI) .CELL  1(1*1  )  DEFINE  THE  LEFT  AND  RIGHT  HANO  EDGES.  * 

C**  RESPECTIVELY.  OF  THE  ITH  CELL*  .tCEJULill)  IS  SET  TQ  * 

C**  -10**75  CELL  1 ( K 1 ♦ 1 )  IS  SET  TO  10**75  BY  THE  PROGRAM  * 

C**  XMAX.XMIN  ARE  MAXIMUM  AND  MINIMUM  X  VALUES*  RESPECTIVELY.  * 

C**  X0AR  IS  THE  COMPUTED  AVERAGE  VALUE  FOR  THE  X  S.  * 

C**  DEV  IS  THE  COMPUTED  STANOARO  DEV  I  AT  ION  FOR.  THE.  X  .S*  _  _ 

C**  Kl  IS  THE  NUMBER  OF  INTERVALS.  * 

CC*  NX  IS  THE  NUMBLR  OF  DATA  POINTS  (sNP)  ♦ 


SUBROUTINE  KAISER  __  .  .  . . 

PURPOSE 

GENERATES  A  SET  OF  KAISER  WEIGHTS. 

CALLING  SEQUENCE  -  -  ~  T 

CALL  KAISER  ( N • SLL • BUF .  IER) 

DESCRIPTION  OF  PARAMETERS  __  _  . 

inputs: 

N  -  NUMBER  UF  SAMPLES 

SLL  -  DESIRED  LEVEL  OF  SIDELOBES  BELOW  MAINLOBE  IN  DB(NEG. 

outputs:  .  . 

HUF  -  ARRAY  FOR  WEIGHTS  * 

IER  -  ERROR  FLAG 

IF  NON-ZERO  •  THEN  INPUT  SLL  WAS  BAO  AND  DEFAULT 
WEIGHTING  IS  CHOSEN. 

COMMON  VARIABLES  USED 
NUN  t_ 

SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

A uthur/prqgrammer 

M  .BURNS 
MAR  1979 
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PROLOGUES  TO  ROUTINES  ON  FILE  1  (CONT. 


) 


cc* 

cc*  note:  a  routine  called  OUT1  (STORED  on  *LCMO.Dl SKL ib*  AS 

CC*  AXOUT1 )  CAN  BE  USED  TO  LIST  THE  RESULTS  GENERATED  BY 

CC*  PRDIST. 

CC* 

CC*  SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  KEUUI RED  1 

CC*  NONE 

CC* 

CC*  AUTHUR/PROGRAMMLR 

CC*  W.  BURNS  . . 

CC*  MAR  19  79 

CC* 

CC** *** ******  *^ ****************** *m** *********************** ***^*^***  m* 


CC*****  ************************  LDRINI  ****************************  *** 

CC* 

CC*  SUBROUTINE  LDRINI 

CC* 

CC*  PURPOSE 

CC*  TO  INITIALISE  FOR  LADDER  PLOTS  BY  DEFINING  PAGE  ORIGIN  AND 

CC*  WRITING  LADDER  LABEL  WITH  TITLE  AND  X.Y  AXIS  DESCRIPTORS. 

CC* 

CC*  CALLING  SEQUENCE 
CC*  CALL  LDRIN! 

CC* 

CC*  COMMON  VARIABLES  USED 

CC*  inputs: 

CC*  /LOR  INF/  XTAG.XUNITS.XLEFT.XRIGHT. XLEN. 

CC*  Y  TAG . YUN ITS .YBOT .Y TOP. YUEN. Tl TLt.RUNGLB 

CC* 

CC*  SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

CC*  PLOT.SYMBUL.CORE  -  CALSUBS 

CC* 

CC*  author/programmer 

CC*  C.M.SSCSESNY 

CC*  FEB  1979 

CC* 

CC*****  ************************************************************** ** 


CC***** ************************  ldrplt  ****************************** 


CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

CC** 


SUBROUTINE  LDRPLT 


PURPOSE 

TO  GENERATE  A  PLOT  OF  Y  VS. 
LADDER  OF  SUCH  PLOTS  * 


X  THAT  WILL  BE  ONE  RUNG  IN  A 


CALLING  SEQUENCE 

call  LDRPLT  ( X, Y • NPTS.RUNGI O ) 


D  E  SCR  I  P  T I  ON 
inputs: 

X 

Y 

NPTS 
RUNG  I 


OF  PARAMETERS 


-  ARRAY  OF  X  DATA  VALUES 

-  ARRAY  OF  Y  DATA  VALUES 

-  NUMBER  OF  POINTS  IN  X.Y 
{>-  IDENTIFIER  FOR  THIS  RUNG 


TO  BE  PLOTTED 
(SUCH  AS  TIME .RECORD 


#.E  TC 


COMMON  VARIABLES  USED 

inputs: 

/LDRXNF/  XLEFT .XR IGHT .XLEN. YBOT  * YT OP .YLEN. RUNGSP. LAB REP 


SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUI RED 
PLOT, NUMBER  -  CALSUBS 
PLT  -  DRAWS  Y  VS.  X 


REMARKS 

usage:  _ 

US!  NG  SUBROUTINE  LDRPLT  REQUIRES  LINKAGE  OF  2  OTHER  SUBROUT  (It 
LDRINI  AND  PLT  (STORED  AS  GNLDRINI  AND  GNPLT  ON  LCMO.  D I  SKLI HJ 
LDRINI  MUST  BE  CALLED  ONCE  BEFORE  1ST  CALL  TU  LDRPLT. 

LDRPLT  IS  CALLED  FOR  EACH  TRACE  THAT  IS  TO  BE  A  RUNG  ON  THE 
l.  ADDER  PLOT.  THE  CALLING  PROGRAM  MUST  T f£N  EXECUTE  A  CALL  TU 
EFPLUT  (ON  CALSUBS)  AFTER  ALL  TRACES  ARE  PLOTTED. 

A  U THOR /PROGRAMMER 
C.M.SZCZESNY 
FEB  1979 

*****************************************************************  *** 
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CC*****  *****4  ******************  PLT  444******  *4** 4  4  4**»****»*4*****  ** 

cc* 

CC*  SUBROUTINE  PLT 

CC* 

CC*  PURPOSE 

CC*  TO  PLOT  THE  I  X.Y  >  COORDINATES  AND  CLIP  THOSE  LINES  THAT  ARE 

CC*  ABOVE  AND  BELOW  THE  PLOT  BOUNDARIES# 

CC* 

CC*  CALLING  SEQUENCE 

CC*  CALL  PLT  (X.Y.NPTS.XMIN • XMAX • XSPAN • YM I N #YM AX# YSPAN > 

CC* 

CC*  DESCRIPTION  OF  PARAMETERS 
CC*  inputs: 

CC*  X  -  ARRAY  OF  X  DATA  VALUES 

CC*  Y  -  ARRAY  OF  Y  DATA  VALUES 

CC*  NPTS  -  NUMBER  UF  POINTS  IN  X#Y  TO  BE  PLOTTED 

CC*  XM1N  *-  MINIMUM  VALUE  FOR  X-AXIS 

CC*  XMAX  -  MAXIMUM  VALUE  FOR  X-AXIS 

CC*  XSPAN  -  LENGTH!  INCHES >  OF  X-AXIS 

CC*  Y  MI N  -  MINIMUM  VALUE  FOR  Y-AXIS 

CC*  Y  MAX  -  MAXIMUM  VALUE  FOR  Y-AXIS 

CC*  Y  SPAN  -  LENGTH!  INCHES )  OF  Y-AXIS 

CC* 

CC*  COMMON  VARIABLES  USED 

CC*  NONE 

CC* 

CC*  SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

CC*  PLOT  -  CALSUBS 

CC* 

CC*  AUTHOR/PROGRAMMER 
CC*  C.M.S^CZESNY 

CC*  FEB  1979 

CC* 

CC*****  4444*4  *************  *  **************  ******  *************  ***  ******  *4> 


CC*****  ****** *4 44  44 4 *4 *********  MAGPLT  **444444444444*4444*4444*444  444 


CC*  4 
CC*  SUBROUTINE  MAGPLT  * 
CC*  4 
CC*  PURPOSE  4 
CC*  TO  GENERATE  A  PRINTER  PLOT  OF  THE  MAGNITUDES  OF  A  2-OIMENSICN4 
CC*  IMAGE  4 
CC*  4 

cc*  calling  sequence  * 

CC*  CALL  MAGPLT  ( V.NR OW » NCOL  •  VSCALE • VTH )  4 
CC*  4 
CC*  DESCRIPTION  OF  PARAMETERS  * 
CC*  INPUTS!  4 
CC*  V  -  ARRAY  CONTAINING  THE  IMAGE  * 
CC*  NRUW  -  1ST  DIMEN.  OF  Y  .(•ELEMENTS  IN  CROSS-RANGE  DIRECJIUNJ 
CC*  NCOL  -  2ND  DIMEN.  OF  V  (  ^ELEMENTS  IN  RANGE  DIRECTION)  * 
CC*  VSCALE-  SCALE  FACTOR  FOR  V  4 
CC*  VTH  -  THRESHOLD  TO  BE  USED  ON  V  * 
CC*  .  ._  .  * 
CC*  COMMON  VARIABLES  USEO  4 
CC*  NONE  4 
CC*  * 
CC*  SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUIRED.-  .  ..  . 

CC*  NONE  4 
CC*  4 
CC*  AUTHOR/PROGRAMMER  4 
CC*  M. BURNS  _  4 
CC*  MAR  1979  4 
CC*  4 


CC***** **************** ************************  4*** 4 ***************** *44 
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PROLOGUES  TO  ROUTINES  ON  FILE  1  (CONT.) 


C.C*****  ************************  OUT  1  ********************************* 

CC*  * 

CC*  ..  SUBR0UUNJl_OU.T1_  _ _ _ 

CC*  * 

CC*  PURPOSE  * 

C**  PRINT  OUTPUT  OF  DATA  STATISTICS  FROM  SUBROUTINE  PRDIST  * 

CC*  _  _ * 

CC*  CALLING  SEQUENCE  ~  * 

CC*  CALL  OUT  I  * 

CC*  * 

CC*  COMMON  VARIABLES  USED  _  _  *  * 

CC*  INPUTSI  '  * 

CC*  /LCPROB/  * 

CC*  NBUCX  -  NO.  OF  OCCURRENCES  FOR  EACH  CELL  * 

CC*  O I  ST  -  CUMMULATIVE  DISTRIBUTIONS  C0RR.LSPQNQI6K*  T.U  J«JGHTL_* 

CC*  HAND  EDGE  OF  EACH  CELL  "  * 

CC*  CELL1  -  DEFINES  THE  EDGES  OF  THE  CELLS  FROM  RIGHT  TO  LEFT* 

CC*  XMAX  -  MAXIMUM  VALUE  * 

CC*  XMIN  -  MINIMUM  VALUE.  __  _  _ 

CC*  XOArt  -  COMPUTED  AVERAGE  '  * 

CC*  DEV  -  COMPUTED  STANDARD  DEVIATION  * 

CC*  K 1  -  NU.  OF  INTERVALS  * 

CC*  NX  -  NO.  OF  DATA  POINTS  * 

cc*  . . . 

CC*  REMARKS  * 

C**  SAMPLE  PHI NT  OUT  COULO  LOOK  LIKE  THIS.  .  * 

C>*  MEAN=  11. 42  SIGMA=  8.6230  M**,  =  JO. 00  M IN=  _ L»OQ  .  NO.  PIS»*=2QJ>-* 

C**  CELL  FREQ.  CUM.  CELL  FREQ .  CUM.  • 

C**  *****  -  0.  0  0.  10.00  -  12.00  20  60.00  * 

C**  0.  -  1.00  O  0.  12.00  -  14.00  10  6b. 00  * 

C**  1.00  -  2.00  20  JO. 00  _l*»Q0  16.00  10 _ JO.OD  * 

C*»  2.00  <-  3.00  10  15.00  16. 00  -  18.00  '  "20  U0.00  • 

C**  3.00  -  4.00  10  20.00  16.00  -  20.00  10  BS.00  * 

C**  4.00  -  5.00  10  25.00  20.00  -  22.00  5  67.50  * 

C*»  5.00  -  6.00  _  10  30.00  _ 22. 00  -  __24jQ0  _  _0  87.50  * 

C**  6.00  -  7.00  10  35.00  24.00  -  26.00  5.  90.00  * 

C»*  7.00  —  8.00  10  40.00  26.00  —  28.00  O  90.0  0  » 

C**  8.00  -  9.00  10  45.00  28.00  -  30.00  0  90.00  * 

.C**  .  _  9.00  10.00  10  50.00  30. 0Q_  **.**** _ 20_IOO.OQ_* 

C»*  A  SKIP  TO  A  NEW  SHEET  PLUS  A  HEADER  LINE  SHOULD  OE  SUPPLIED  * 

C**  BY  THE  USER  IN  THE  MAIN  CALLING  PROGRAM.  4 

CC«  _  _  3 

CC*  SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUI  RfcD  -  •  -  , 

CC*  NONE  » 

CC*  4 

CC*  AUTHOR /PROGRAMMER  .  _  __  _  _  _  _  _* 

CC*  BURNS/SZCZESNY  "  """  -  4 

CC*  MAR  1979  4 

CC*  » 

CC****»  ********  **********************  ****♦*.«*  ****.*********  ***44414.4*  *4 
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cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 


******** ******************************** ********** **************** 
**  IMAG2  •* 

X-BANO  CONTOUR  PLOT  PROGRAM 

PURPOSE 

TO  READ  A  TAPEL  OF  WIDEBAND  RADAR  DATA*  CONVERT 
THESE  DATA  INTO  MAGNf TUDE  VALUES  AND  PRODUCE  TWO 
DIMENSIONAL  CONTOUR  IMAGES  OF  THE  CONVERTED  DATA 

INPUTS 

cards: 

VARIABLES  TO  CONTROL  RE AD I NG .PR  INTI NG •  PLOTTING 

method: 

3  namelist  statements  and  1  formatted  reao  statement  are 

USED  TO  INPUT  PARAMETERS  AS  SHOWN  BELOW: 

NAMELIST  /DO INOX/  -  PARAMETERS  TO  CONTROL  READING  OF 
THE  INPUT  TAPE 

NREC1  -  THE  RECORD  NUMBER  (POSITION  >  OF  THE  INITIAL 
RECORD  TO  BE  READ  FOR  EACH  IMAGE  * 

NREC2  “  THE  RECORD  NUMBER  (POSITION)  OF  THE  LAST 
RECORD  TO  BE  READ  FOR  EACH  IMAGE 
N 1 MAGS  -  THE  TOTAL  NUMBER  OF  IMAGES  ON  THE  TAPE 
NIMAGP  -  THE  NUMBER  OF  THESE  IMAGES  TO  BE  PROCESSED 
DURING  THIS  RUN 

FORMATTED  READ  (1000)  -  PARAMETER  TO  CONTROL  WHICH  IMAGES 


CC* 

CC* 

CC* 

CC* 

CC* 

CC* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 


ARE  PLOTTED*  USED  IN  CONJUNCTION 
WITH  THE  NAMELIST  DOINDX 

IMAGNO  -  AN  ARRAY  OF  IMAGE  NUMBERS*  IN  ASCENDING  ORDER* 
OF  THOSE  IMAGES  TO  BE  PLOTTED  FOR  A  GIVEN  TAPE* 
SEE  COMMENTS  BELOW  IN  THE  SOURCE  FOR  MORE  INFO 


NAMELIST 
I  PR  I  NT 
IF  0 
IF  I 
IF  2 
I  PLOT 
IF  0 
IF  1 
ISOALE 
IF  0 


IF 


NLEV 

CTHV 

CINCR 

SIZEX 

SIZEY 

XUNITS 

YUNITS 

XM1N 

YMIN 

yorign 

XOR1GN 

LABELX 

LABELY 

LSLSIZ 


/1FLAGS/  -  PRINT  AND  PLOT  FLAGS 

-  ASSUMES  A  VALUE  OF  0 •  1  OR  2 

-  NO  OUTPUT  TO  PRINTER  OF  MAGNITUDE  OR  PHASE  VALUES’ 

-  MAGNITUDE  VALUES  ARE  OUTPUTTED  TO  PRINTER 

-  MAGNITUDE  AND  PHASE  VALUES  ARE  OUTPUTTED 

-  ASSUMES  A  VALUE  OF  0  OR  I 

-  NO  CONTOUR  PLOTTING  OF  MAGNITUDE  VALUES  IS  DONE 

-  CONTOUR  PLOT  THE  MAGNITUDE  VALUES 

-  ASSUMES  A  VALUE  OF  0  OR  I 

-  USER  SUPPLIED  INPUTS  FOR  SIZE  OF  X  AND  Y  AXIS 
ISIZEX.SIZEY)  AND  NUMBER  OF  UNITS  OF  X  AND  Y  PER 
INCH  FOR  EACH  AX 1 S( XUNI TS* YUNI TS )  ARE  PASSED  TO 
THE  PLOT  SUBROUTINE  UNCHANGED* 

-  USER  SUPPLIED  INPUTS  FOR  SI ZEX .XUNI TS*YUNI TS 
ARE  IGNORED  AND  THE  PROGRAM  CALCULATES  VALUES 
FOR  THE  3  VARIABLES  BASED  ON  THE  RATIO  OF  CROSS 
RANGE  CELLS  TO  DOWN  RANGE  CELLS*  SEE  SOURCE 
COMMENTS  FOR  MORE  INFORMATION 

NAMELIST  /IMPLOT/  “  PARAMETERS  TO  CONTROL  THE  CONTOUR  PLOT 
PROCEDURE* 

NUMBER  OF  CONTOUR  LEVELS  DES1REO 

CONTOUR  THRESHOLD  VALUE! SMALLEST  VALUE  CONTOURED) 
INCREMENT  BETWEEN  CONTOUR  LINES 
SIZE  OF  X  AXIS  IN  INCHES 
SIZE  OF  Y  AXIS  IN  INCHES 

NUMBER  OF  UNITS  OF  X  PER  1NCH!E*G*  1.SM/IN) 

NUMBER  OF  UNITS  OF  Y  PER  INCH 
INITIAL  VALUE  FOR  X  AXIS 
INITIAL  VALUE  FOR  V  AXIS 
POSITION  OF  X  ORIGIN  IN  INCHES 
POSITION  OF  Y  ORIGIN  IN  INCHES 
ARRAY  OF  CHARACTERS  FOR  X  LABEL 
ARRAY  OF* CHARACTERS  FOR  Y  LABEL 

NUMBER  OF  CHARACTERS  IN  X  OR  Y  LABEL! SAME  SIZE) 


TAPE  • 

RAW  WIOEBANO  PULSE  DATA  TO  BE  CONVERTED  TO 
MAGNITUDE  VALUES  FOR  PLOTTING 
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cc* 

CC*  OUTPUTS 

CC*  STSPR1NT  -  HEADER  record  for  each  image,  also  phase 

CC*  AND  MAGNITUDE  VALUES  IF  PRINT  FLAG  IS  "ON** 

CC*  PLOTTER  -  MAGNITUDE  VALUES  ASSOCIATED  WITH  EACH  IMAGE 
CC*  ARE  CONTOUR  PLOTTED  ALONG  WITH  APPROPRIATE 

CC*  HEADER  INFO  USING  THE  VERSATEC  PLOTTER 

CC*  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
CC* 

cc*  internaliibm  supplied)  fortran  functions  _ 


cc*  A I  MAG  -  EXTRACTS  IMAGINARY  PORTION  OF  COMPLEX  VAR  •* 

CC*  REAL  -  EXTRACTS  REAL  PORTION  OF  COMPLEX  VARIABLE  •< 

CC*  ATAN2  -  TANGENT  FUNCTION  .RETURNS  A  A  QUADRANT  RESULT*. 

CC*  CABS  -  ABSOLUTE  VALUE  OF  COMPLEX  VARIABLE  *< 

CC*  INT  -  1NTERGER  PORTION  OF  FLOATING  POINT  VARIABLE  •< 

CC*  ALOGIO  -  LOG  FUNCTION  * 

CC*  • 

CC*  EXTERNAL  SUBROUTINES  • 

cc*  MCCNTR  * 

CC*  PURPOSES  PROOUCE  CONTOUR  PLOT  OF  DATA  • 

CC*  IN  A  2-DIMENSIONAL  ARRAY  * 

CC*  REQUIRES:  PLOT. AXIS. SYMBOL. NUMBER  • 

CC*  ARGUMENTS  ( 2-0  I MENS  I ON AL  ARRAY .PLOT  • 

cc*  parameters )-see  description  of  sue  • 

CC*  MCCNTR  FOR  MORE  DETAILS  * 

CC*  • 

CC*  authorsprogrammer  • 

CC*  CAL  SPAN  •' 

cc*  feb  1979  * 

CC*  P.M  .MCMAHON  *• 

CC*«*****«* ******** ********* *********************9*****^****************' 
cc*********** ******************  MCCNTR  a********************************' 
CC** 

CC**  SUBROUTINE  MCCNTR 

cc**  '7 

CC**  PURPOSE  *' 

CC**  TO  GENERATE  PLOTTED  CONTOUR  IMAGES  OF  A  2-0 IMENSIONAL  ARRAY  *' 

CC**  OF  DATA  POINTS 

cc** 

CC**  CALLING  SEQUENCE  * 

CC**  CALL  MCCNTRIA.NI.NX.NV.C.NC.XMIN.DX.ALX.TNIN.DY.ALT.XI.VI.  * 

CC**  xlab.ylab.nl. I PXYt  * 

cc**  • 

CC**  DESCRIPTION  OF  PARAMETERS  • 

CC**  inputs:  * 

Cc**  a  -  2  dimensional  array  of  points  to  be  plotted  • 

CC**  N1  -  NUMBER  OF  ROMS  IN  THE  ARRAY  IN  THE  CALLING  PROGRAM  • 

CC**  NX  -  NUMBER  OF  X  VALUES  • 

CC**  NY  -  NUMBER  OF  Y  VALUES  • 

CC**  C  -  I  DIMENSIONAL  ARAY  OF  CONTOUR  LEVELS  TO  BE  PLOTTED  * 

CC**  NC  -  NUMBER  OF  CONTOUR  LEVELS  TO  BE  PLOTTED  * 

CC**  XMIN  -  X  VALUE  OF  FIRST  POINT  IN  THE  ARRAY  I  USED  AS  • 

CC**  THE  INITIAL  VALUE  OF  THE  X-AXIS)  * 

CC**  OX  -  NUMBER  OF  UNITS  OF  X/ INCH  IN  THE  PLOT  (USED  AS  THE  * 

cc**  interval  for  the  x  axis)  * 

CC**  ALX  -  LENGTH  OF  X  AXIS  IN  INCHES  * 

CC**  TMIN  -  Y  VALUE  OF  FIRST  POINT  IN  ARRAY  A  (USED  AS  THE  • 

cc**  initial  value  of  the  v-axisi  • 

CC**  DY  -  NUMBER  OF  UNITS  OF  X/1NCH  IN  THE  PLOT  (USED  AS  THE  • 

CC**  INTERVAL  FOR  THE  Y-AXIS)  * 

CC**  ALV  -  LENGTH  OF  Y-AXIS  IN  INCHES  * 

CC**  XI  -  COORDINATE  IN  INCHES  OF  THE  X  ORIGIN  (ORIGIN  RESET  * 

CC**  FROM  0.0  TO  XI. VI)  * 

CC**  Y 1  -  COORDINATE  IN  INCHES  OF  THE  T  ORIGIN  * 

CC**  XLAB  -  ARRAY  OF  CHARACTERS  TO  BE  USED  AS  THE  X  LABEL  * 

CC**  VLAB  -  ARRAY  OF  CHARACTERS  TO  BE  USED  AS  THE  Y  LABEL  * 

CC**  NL  -  NUMBER  OF  CHARACTERS  IN  EITHER  LABEL  (NOTE:  LENGTH  * 

CC**  OF  XLAB  «  LENGTH  OF  VLAB)  * 

CC**  IPXY  -  PLOT  FLAG.  ALMAYS  SET  AT  2  FOR  THIS  TYPE  OF  DATA  .  • 

CC**  IPXY  a  I-  PLOT  IS  SVMETRIC  ABOUT  X-AX  IS.  CORRESPOND  I  N«  • 

6C**  POINTS  PLOTTED  IN  UPPER  QUADRANTS  * 

CC**  IPXY  ■  2-  PLOT  ONLY  THOSE  LEVELS  FROM  THE  GIVEN  * 

CC**  POINTS — PLOT  THE  DATA  AS  GIVEN  • 


cc** 

inputs: 

cc** 

A 

— 

cc** 

N1 

— 

cc** 

NX 

— 

cc** 

nv 

— 

cc** 

c 

— 

cc** 

NC 

— 

cc** 

XMIN 

— 

cc** 

cc** 

ox 

cc** 

cc** 

ALX 

_ 

cc** 

YMIN 

— 

cc** 

cc** 

DV 

_ 

cc** 

cc** 

ACT 

cc»* 

XI 

• 

cc** 

cc** 

Y1 

cc** 

XLAB 

— 

cc** 

VLAB 

— 

cc** 

NL 

— 

cc** 

cc** 

IPXY 

_ 

cc** 

cc** 

cc** 

cc** 

cc** 

cc** 

cc** 

REMARKS 

THE  RC<X 

THE  REQUESTED  CONTOUR  LEVELS  ARE  INTERPOLATED  FROM  THE  DATA 
POINTS  AM)  THEN  MRITTEN  TO  A  DATA  SET  AND  SORTED  BEFORE 
BEING  PLOTTED  ON  THE  VERSATEC 


CC**  POINTS  AM)  THf 

CC**  BEING  PLOTTED 

CC** 

CC**  AUTHOR/PROGRAMMER 
CC**  E. BECKER 

CC**  FEB  1979 
CC«*  P.M. MCMAHON 


>••••••*••••  •*••*****•*•**•*••**•**• 
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PROLOGUES  TO  ROUTINES  ON  FILE  3 


CC ***********  ******  ************  NO  AND  ****♦********♦*♦***♦♦**♦*♦*♦**111 

CC* 

CC*  AXD  NARROW  BAND  PLOT  PROGRAM 

CC*  ' 

CC*  PUR  post 

CC*  ro  READ  SELECTED  RECORDS  UK  A  TAPE  OK  NARROW  BAND  RADAR  DATA 

CC*  WHERE  EACH  RECORD  IS  COMPOSED  OF  RADAR  CROSS  SEC T I ON , T I  ME  • 

CC*  ASPECT  ANGLE  AND  POKJECTfcO  ASPECT  ANGLE  AND  TMtN  PERFORM  THE 

CC*  FOLLOWING  OPERATIONS 

CC*  1)  SMOllT  H  THE  DATA  ACCORDING  TO  USER  SUPPLIED  PARAMETERS 

CC*  2)  CALCULATE  DESCRIPTIVE  STATISTICS  ON  THE  RCS  DATA 

CC*  3)  PLOT  RCS  VS  TIME  OR  RC5  VS  ASPECT  ANGLE  Uk  RCS  VS  *  "" 

CC*  PROJECTED  ASPECT  ANoLL 

CC* 

CC*  DESCRIPTION  UF  PARAMETERS 
CC*  INPUTS 

cc*  cards: 

CC*  VARIABLES  TO  CONTROL  SMOOTHING.PLOTTlNG.STAT ISTICS. TAPE  READING 

CC*  METHOD- 

CC*  4  NAMELIST  READ*  STATEMENTS  ARE  USED  TO  INPUT  PARAMETERS 

cc*  as  shown  oelow: 

CC* 

Cc*  NAMELIST  /RtADTP/-  PARAMETERS  TO  CONTROL  READING  OF  THE 

CC*  INPUT  TAPE 

CC*  NRECl  THE  RtCORU  NUMBER  UK  THE  INITIAL  RECORD  TO  BE  READ 

CC*  FROM  THE  TAPE 

CC*  NRECa  -  THE  RECORD  NUMBER  OF  THE  LAST  RECORD  TO  BE  READ 

CC*  FROM  THE  TAPc 

cc* 

CC*  NAMELIST  /IFLAGS/  -  CONVERS I ON . ST  A T l S T l CS . PLO T  AND  TYPE 

CC*  »  UF  KCS  UNITS  FLAGS 

CC*  ICNVRT  -  ASSUMES  A  VALUE  OF  0  OR  1 

CC*  IF  0  -  NO  CONVERSION  OF  RCS  DATA  BETWEEN  SM  AND  OBSM 

CC*  IF  1  -  CALL  SUBROUTINE  TO  CONVERT  RCS  ARRAY  FROM 

CC*  RCS/SM  TO  RCS/OUSM 


CC* 

CC* 

CC* 

cc* 

cc* 

CC* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 


I  STATS 
IF  o 

IF  1 
I  plot 
IF  0 
IF  1 

IF  2 
IF  3 

T  YPRCS 

IF  o 

IF  1 


ASSUMES  A  VALUL  OK  0  OR  1 

DU  NOT  CALCULATE  DESCRIPTIVE  STATISTICS  ON  The  " 
RCS  DATA  ARRAY 

CALL  SUBROUTINE  TO  CALCULATE  Dt SCRCR 1 PT I VE  STATS 
ASSUMES  A  VALUE  OF  0,1.2  OR  3 
■  NO  PLOTTING  OF  ANY  DATA  IS  PERFORMED 
PRODUCE  A  PLOT  OF  T | ME < X )  VS  KCs! Y )  VALUES 
PROOUCE  A  PLUT  OF  ASPECT  ANGLE! X)  VS  RCS  VALUES 
PRODUCE  A  PLUT  UK  PROJECTED  ASPECT  ANGLE! X>  VS 
RCS(Y)  VALUES  • 

ASSUMES  A  VALUE  OK  0  OR  1 

RCS  IN  SM  UNITS,  THEREFORE  NO  AUTOMATIC  SYMETRICAL 
SCALING  OF  RCS! Y )  AXIS  ABOUT  0 

RCS  IN  OUSM  UNITS.  THEREFORE  AUTOMATICALLY  SCALE 
SYME  f R I CALLY  ABOUT  0  RCS 


NAMELIST  /I PULSE/  -  PARAMETERS  TO  CONTROL  SMOOTHING  AND 

'data  output 

NPLSES  -  THE  NUMBER  UF  PULSES  TO  BE  AVERAGED  TO  PRODUCE 
ONE  RCS  DATA  POINT  FOR  PLOTTING  OR  STATISTICS 
NTHUUT  -  CONTROLS  HOW  MANY  SMOOTHED  DATA  POINTS  ARE 
OUTPUTED  TO  PRINTER  FOR  CHECKING  PURPOSED 


NAMELIST  /STATS/  -PARAMETERS  TO  CONTROL  THE  FREQUENCY 
DISTRIBUTION  THAT  IS  UUTPUTCD  BY 
"*"■  "  '  THE  STATISTICS  ROUTINE 

tape: 

NARROW  BAND  RADAR  RECORDS  COMPOSED  OF  RADAR  CROSS  SECTION. 

TIME  OF  GOSFRVaTIuN  AND  Two  ASPECT  ANGLES 
outputs: 

SY SPRINT  -  HEADER  RECORD  FOR  THL  F I LE *  A  RECORD  OF  DATA  FROM 
THE  TAPE  AT  A  USER  SPECIFIED  INTERVAL.  AND 
~  DESCRIPTIVE  STATISTICS  INC.  FftEuUENCY  01 STR IBUT ION* IF 

•  THE  I  STATS  FLAG  IS  "ON" 

PLOTTER  -  A  LINE  PLOT  OF  RCS(DUSM  OR  SMI  VS  TIME. ASPECT  ANGLE  OR 
PROJECTED  ASPECT  ANGLE  IF  I  PLOT  =  1.2  OR  3  RESPECT  I VEL 
USING  CALSPANS  LAlCOMP  PLOTTER 
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PROLOGUES  TO  ROUTINES  ON  FILE  3  (CONT.) 


CC*  SUURUUT  INt  aN:i  FUNCTION  SUBPROGRAMS  REUVMREO  > 

CC*  SMOOTH  -  PERFORMS  SMOOTH  1.40  UPt^AT  iiH  ON  THE  PCS  DATA 

CC*  CONVRT  -  CONVERTS  AN  ARRAY  OF  RCS/SM  DATA  TO  RCS/uUSM  OAT  A  '  1 

CC*  PROIST  -  CALCULATE  MEAN  •  ST  AND  AMO  DL V I  AT I  ON  »MIN*MAX  A\ll|f 

CC*  rwEUUENCY  1)1  STR  I  OUT  ION  OF  AN  ARkAY  OF  KCS  DATA  ‘ 

cc*  uori  -  uurpurs  results  calculated  in  prdist  > 

CC*  PlTTME  -  DMA  /S  A  PlIT  Or  RCS  VS  THE  UN  THE  CALCUMP 

CC*  PL.  T  A  SP  *  DPA.V.J  PLOT  OF  KC5  VS  ASPECT  ANGLE  OR  RCS  VS  PROJtC  TED  1 

CC*  ASPLwT  ANGLr,  ON  CALOJMP 

CC* 

CC*  AUTHOR^PROGRAMMER  1 

CC*  P.M.MCMAHQN 

CC*  APR  1V79 

CC*  C  AC  SPAN  CORPORATION 

CC* ********  ********  *********** ***********  ************** **************** 

cc*  _ 

cc* **********  ******************  ^Miiu r m  ****************************** 

cc* 

CC*  SUBROUTINE  SMOOTH 

CC* 

CC*  PURPOSE 

cc*  to  h«.au  opuses  records  of  data  from  the  narrow  band  tape  •calculate 

CC*  THE  MEmN  OF  r  Hr.  R  C'J  DATA  Pul  NTs  AAi)  SELECT  T>1r.  MEDIAN  TIME*  ASPECT 

CC*  ANCLi-  AND  PROJECTED  ASPECT  VAlUlS 

CC* 

cc*  c  all  I  fKi  Seoul  Net  _ . _ 

CC*  CALc  SMOOTMlNPUSEStRCSVAL,  TMEVAL*ASPCTl*ASPCr*.»,cOF) 

CC* 

CC*  DESCRIPTION  of  PARAMETERS 

cc*  iNPury 

cc*  argo/mlnt  list: 

CC*  NPUSeS  -  USER  DEFINED  VALUE  FOR  the  NUMUlK  OF  RECORDS  TO  BE 

CC*  Rt  AU  ANU  A  YE  RACED 

cc*  tape: 

CC*  NARROW  HAND  RADAR  TAPE  *ITH  RECORDS  COMPOSED  OF  A  RCS  VALUE* 

CC*  TIME  VALUE  AMU  d  ASPECT  ANCLES 

CC*  OuTpuTs 

CC*  ARGUMENT  LIST: 

CC*  RCSVAL  -  ME  A  \|  OF'  THE  RCS  VALUtiSl  AFTER  CONVERTING  UNITS  FROM  DJSM 

cc*  ro  sm)  for  nplsls  records 

CC*  TMtVAL  -  MEDIAN  TIME  VALUE  OF  THE  NPLSES  RECORDS 

CC*  ASPCTI  -  MEDIAN  ASPECT  ANCLE  OF  THE  NPLSES  RECORDS 

CC*  ASPCTii  -  MEDIAN  PROJECTED  ASPECT  ANGlF.  Of  THE  NPLSES  RECORDS 

CC*  LOF  -  END  OF  FILE  FLAG*  SET  TO  I  IF  EOF.  IS  RbACMeD  WrllLe 

CC*  READING  THE  TAPE  WIThIN  SMOOTH 

CC* 

CC*  Rt MARKS 

CC*  THE  MAIN  PROGRAM  ADVANCES  THE “TAPE  TO  THE  STARTING  RECORD  AND* 

CC*  TESTS  FOR  THE  END-UF-FIlE  CONDITION.  IF  EOF-1 •  THE  OTHER  VALUES 

CC*  BEING  RETURNED  ARE  IGNORED 
CC* 

CC*  subroutine  AND  FUNCTION  SUBPROGRAMS  REUUIRED 

CC*  NONt 

CC* 

CC*  AUTHOR/PROGRAMMER 

CC*  P.M. MCMAHON 

CC*  APR  1R79 

CC*  CAL SPAN  CORPORATION 

CC* ******** **  ********************************************************** 
CC* ******************* *********  CONVRT  *******Mt****«*«***4*u»m*« 

CC* 

CC*  SUBROUTINE  CONVRT 

CC* 

CC*  PURPOSE 

CC*  CONVERTS  AN  ARRAY  OF  KCS(SM)  DA TA  TO  HCS(DBSM) 

CC* 

CC*  CALLING  SEQUENCE 
CC*  CALL  CON  VRT IN*  RCSQOB  ) 

CC* 

CC*  DESCl?  IP T  I  UN  OF  PARAMETERS 

CC*  INPUTS 

CC*  N  -  NUMBER  OF  DATA  POINTS  IN  THE  RCS  ARRAY ( MAX= 1 0000 > 

CC*  RCSQD0  -  THE  RCS  ARRAY  IN  OR»I  G I  NAL  (  SQ  •  M)  UNITS 

CC*  OUTPUTS 

CC*  RCSQOB  -  THE  RCS  ARRAY  IN  C ON VKR TED I DBSM >  UNITS  Fqr  PLOTTING 
CC* 

CC*  REMARKS 

CC*  THF.  ORIGINAL  UNITS  ON  THt  TAPE  IS  OBSM  HOWEVER*  THE  SMOOTH  SUB 

CC*  CONVERTS  TO  SM  BEFORE  AVERAGING*  THIS  ROUTINE  CONVERTS  BACK  TO 

CC*  DQSM  FOR  PLOTTING  OR  STATISTICS 


-HIS  r&az  is  waiw&Bijr 

i'hvu  jojj  raojr^rjL)  rt.  iy>,q  ^  .. 
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cc* 

cc* 

cc* 

cc* 

cc* 

cc* 

cc* 


suurout Inc  and  function  subprograms  required 
1 DM  SUPPLIED  tOG  FUNCTION  ALOClO 


AUTHOR/ PROGRAMMER 
P.M  •MCMAHON 
APR  1979 

CAL  SPAN  CORPORATION 

CC*  ******  ****  ***********************  A**  *4>  ***************************** 

C******************************************^*******  *************  *+***« 

c*  SUBROUTINE  PRO  I  ST 
C* 


c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

t: 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c***< 


PURPOSE  '  '  — 

liVALUAT  I  UN  uP  DATA  STATISTIC^ 

1 )  MEAN  VALUE 

2)  STANDARD  DEVIATION 

3)  MAXtMUA 
A)  MINIMUM 

b)  Hi  STOGRAMCPRCOUENC  l£S  UF  OCCURRENCE) 
b)  CUMULATIVE  DISTRIBUTION 

USAGE 

CALL  PROISTIX  .NP.A.B.C.O.t ) 

WHERE 

X  IS  THE  NAME  FOR  THE  DATA  ARRAY 

np  is  the  number  of  data  points 

DEFINE  CELL  bIZcS  AS 

A(B)C(D>E  -  (  THAT  IS.  A  TO  C  IN  INCREMENTS  OF  SIZE  B. 

then  ro  E  IN  INCREMENTS  OF  SIZE  D). 

FOR  A ( B ) C •  SET  £=D=C  UR  SET  E=  0.0. 

REMARKS 

RESULTS  ARE  LEFT  IN  LABELED  COMMON  LCPRUD.  DEFINLD  BY 
COMMUN/LCPRUli/  NBUCKt  lOl).DlST(lGl)  .CELL  i  <  I  02  )  .X MAX, AMIN. 
XBAR.DLV .Kl .NX 

NBUCKI  I)  IS  THE  Nd.  OF  OCCURRENCES! I NTEGER I  FOR  THE  ITH  CcLl 
DISTU)  4S  THE  CUMULATIVE  D  I  S  TR  I  BUT  ION  (  Fl.uA  T  I  NO  PUlNT  VALUE 
INTERVAL  0.0  TO  1.0)  CORRESPOND  ING  ft)  THE  RIGHT  HAND 
EOGE  (JF  THE  ITH  CELL. 

CELL  I <  I )  .CelL 1(1*1)  DEFINE  THE  LEFT  AND  RIGHT  HAND  EDGES, 
RESPECTIVELY.  OF  THE  ITH  CELL.  ( CELL  1(1)  IS  SET  10 
-10A*7S  CELL l (K1+ 1 )  IS  SET  TO  1 0**7S  BY  THE  PROGRAM 
XMAX.XMIN  ARE  MAXIMUM  AND  MINIMUM  X  VALUES.  RESPECTIVELY. 
XBAR  IS  THE  COMPUTED  AVERAGE  VAi.Ub  FOR  THl  x  S. 

DEV  IS  THE  COMPETED  STANDARD  DEVIATION  FOR  THE  X  S. 

XI  IS  THE  NUMBER  OF  INTERVALS. 

nx  is  the  number  of  data  points  <  =  np) 

SUBROUTINE  AVAILABLE 

THIS  ROUTINE  IS  AN  OPTIONAL  OUTPUT  RUUTINh  COMPATIBLE  WITH 
THE  PRECEDING  PRDIST  ROUTINE. 

CALLING  SEUUENCE  -  CALL  OUT  I 

A  SKIP  TO  A  NEW  SHEET  PLUS  A  HEADER  LINE  SHOULD  HE  SUPPLIED 
THE  USER  IN  THE  MAIN  PROGRAM. 


*' 

*■ 

• 

IN* 


********  ************4<**4 


BY* 
* 
* 
♦♦A 


C  ******  ******  *******41************************************************** 


C» 

sun ROUT INE 

OUT  1 

C* 

PURPOSE 

C* 

PRINT 

OUTPUT  OF  DATA 

statistics 

FRUM  SUBRUUTINE 

PRD  I  ST 

CA 

USAGE 

C* 

CALL  OUT ] 

c* 

REMARKS 

c* 

SAMPLE 

PRINT  OUT  COUlD  LOOK 

LIKE  THIS. 

c* 

MEAN  =  11.42  S I  (»MA=  8.b23b  MAX 

»  = 

30.00  MIN 

=  1.00 

NO.  MTS. =200 

c* 

CELL 

FREU. 

cum. 

CELL 

FREQ. 

CUM. 

c* 

A  A*  AA 

0. 

O 

0. 

10.00 

— 

12.00 

20 

O0.00 

c* 

0. 

1  .  UO 

0 

0. 

1  2  •  oO 

— 

14.00 

1  0 

ob, 00 

c* 

1. 00 

2.00 

20 

10.00 

14.00 

— 

lb.  00 

1  U 

70.00 

c* 

2.00 

3.00 

10 

1  S  •  0  0 

lb. 00 

— 

1  8  •  0  0 

20 

80.00 

c* 

3.  00 

4.00 

1  0 

20.00 

1H.00 

— 

20.00 

10 

OB.OO 

c* 

4.  00 

b.  00 

lu 

2S-00 

20.00 

— 

22.00 

S 

tt  /.SO 

C  A 

b.  00 

-  0. 00 

10 

30.00 

22.00 

— 

24.00 

0 

87 ,50 

c* 

6.  00 

7.00 

I  0 

35.  Oo 

24.00 

— 

2b.  00 

b. 

90,00 

CA 

/•OO 

<3.00 

1  0 

40. 00 

2b  iOu 

- 

28.00 

0 

RO  .00 

CA 

a.  oo 

9*0  0 

1  0 

4S.  0(1 

2  a.  oo 

— 

Jo.  00 

o 

90  .00 

CA 

9.  OO 

10.00 

10 

SO. 00 

JO.  00 

— 

****** 

20 

1  00.00 

CA 

CA 

A  SKIP 

TO  A  NEW 

SHEeT 

PLUS  A 

HEADLR  LINE 

SHOULD 

BE  SUPPLIED 

CA 

BY  THE 

USER  IN 

the  main  calling 

PROGRAM 

. 

C* ***********  ************  *********************************************** 


THIS  FAGS  IS  FF.~T  V’-M-TTY  FSAC1 
f/wV.  COFi  bl>G 


ICABI 
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CC***** ****** ************ ******  PLTTMt  I,****************************** 

cc*  * 

CC*  SUBROUTINE  PL  TIME  * 

CC*  *'  . .  '■  “  ■  . . .  * 

CC*  PURPOSE  * 

CC*  DRAW  A  LINE  PLOT  OF.  RADAR  CROSS  SECT  ION<  DBSM  OR  SM  >  VS  TIME  FOR  * 

CC*  REAL  OR  DERlVEDtFROM  WIDEBAND)  NARROW  BAND  RADAR  DATA  * 

cc*  * 

CC*  CALLING  SEQUENCE  * 

CC*  CALL  PLTTMEI  T I Mt • RCSD T A , NDPT S .NREC 1 . NKEC2 .NPULSE . NOB J .FLG >  * 

CC*  * 

CC*  DESCRIPTION  OF*  PARAMETERS’  *• 

CC*  INPUTS  *< 

CC*  ARGUMENT  LIST;  *< 

CC*  TIME  -  ARRAY  OF  TIME  DATA  POINTS  TO  BE  PLOTTED  ** 

CC*  RCSOTA  -  ARRAY  *OF'RCS  DATA  POINTS  TO  BE  PLOTTED  ~  *t 

CC*  NDPTS  -  ARRAY  OF  RCS  DATA  POINTS  IN  TIME  UR  RCSOTA  *< 

CC*  NR  EC  1  -  POSITION  OF  FIRST  RECURD  READ  FROM  TAPE  -  TO  PROVIDE  *• 

CC*  IDENTIFICATION  INFORMATION  ON  THE  PLOTTED  OUTPUT  *. 

CC*  NREC2  ‘=r  POSITION' "OF  CAST  RECORD  READ  FROM  TAPE  -  TO  PR0V1UE - 

CC*  IDENTIFICATION  INFORMATION  AS  ABOVE  * 

CC*  NPULSES  -  SMOOTHING  INTERVAL  SELECTED  BY  USER  -  FOR  PLOT  * 

CC*  I DENT IF  I  CAT  I ON  * 

cc*  nohj  -  object  Number  -  for  plot  identification  . ~* 

cc*  flg  -  .a  flag  that  controls  automatic  symetrical  scaling  of  * 

CC*  RCS  DATA  IY  AXIS)  ABOUT  01. E.  IF  F|_G=1.  SCALING  IS  * 

CC*  AUTOMATIC,  if  not  |.  UStR  MUST  SUPPLY  BEGINNING  AND  * 

CC*  ENDTNG  POINTS  OF  THE  Y  AX  1 5  * 

cc*  outputs:  * 

CC*  LINE  PLOT  OF  RADAR  CROSS  SECTION  VERSUS  TIME  WITH  APPROPRIATE  * 

CC*  IDENTIFICATION  INFORMATION  ON  CALSPANS  CALCOMP  PLOTTER  * 

CC*  .  * 

CC*  SUBROUTINE  AND  FUNC T I  UN  SUBPROGRAMS  REQUIRED  * 

CC*  CALSPAN  PLOT  SUURUUT INES  * 

CC*  PLOT  -  BASIC  PLOT  SUB  USEO  BY  HIGHER  LEVEL  PLOT  ROUTINES  ,* 

CC*  AXIS  -  DRAWS  X  OR  "Y  "ASCI  S  WITH  TIC  MARKS  AND  LABELS  ? 

CC*  SYMBOL  -  DRAWS  Ai-PHA  CHARACTERS  OF  SPECIFIC  SIZE  AND  LOCATION  • 

CC*  NUMBER  -  DRAWS  NUMERIC  INFORMATION  OF  SPECIFIED  SIZE  AT  • 

CC*  A  SPECIFIED  LOCATION  • 

CC*  LINE  -  DRAWS  A  STRAIGHT  LTNE  BETWEEN  POINTS  DEFINED  UY” . — - - H 

CC*  SUCCESSIVE  SLEMENTS  OF  X  AND  Y  ARRAYS  *t 

CC*  IBM  SUPPLIED  FORTRAN  FUNCTIONS  *t 

CC*  INT  -  RETURNS  INTEGER  PUR T ION  OF  FLOATING  PLDlNT  NUMBER _ *l 


"RETURNS  MAXIMUM  VALUE  OF  ARC  OF  FLOATING  PT  NUMBERS 
"RETURNS  “MINIMUM  VALUE  OF  ARG  OF  FLOAT TKIG  "PT  'ROMBERG"' 


CC*  AMAXI 

CC*  AM  INT 

CC* 

CC*  AUTHOR/ PROG  RAMMER  •< 

CC*  P .M .MCMAHON 

CC*  APR  1 9T<? 

cc*********************************************************^************ 


..  ,.rt  l  V*A,“ 
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TABLE  B3 


PROLOGUES  TO  ROUTINES  ON  FILE  3  (CONT.) 
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4  SUBROUTINE  PL  T ASP 

PURPOSE 

DRAW  A  LINE  PLOT  OK  RADAR  CROSS  SEC  T  ION  (OOSM  OR  SM  )  VS  ASPECT 
ANGLL  UR  PROJECTED  ASPECT  ANGLt  FOR  SELECTED  PORTIONS  OF  NARROW 
BAND  RADAR  TAPE 

CALLING  SEQUENCE 

CALL  PLT  ASP  (  ANGLE  •  RC  SDT  A  .  NDP  TS  t  NREC  1  •  NRE  C£  •  NPULSE  *NOB J . FLG  I 
REMARKS 

IN  ORDER  FUR  THIS  PLOT  SUBROUTINE  TU  WuRK  PROPERLY  THt  DATA  THAT 

IS  PAS  SEO  TO  IT  MUST  t3E  IN  ASCENDING  OR  DESCENDING  ORDER 

THEREFORE  UTMOST  CARE  MUST  BE  USED  IN  SPECIFYING  NREC I  AND  NREC 
TO  INSURE  THAT  THE  ANGLE  CHOSEN  IS  INCREASING  OR  DECREASING 
In  one  direction  only 

DESCRIPTION  OF  PARAMETERS 
INPUTS 

ARGUMENT  LIST: 

ANGLE  -  ARRAY  OF  ASPECT  AnGuE  OR  PROJEC TED  ASPECT  ANGLE  DATA 
POINTS  TO  TJE  PLOTTED 

RCSOTA  -  ARRAY (OF  RCS  DATA  POINTS  TO  Dt  PLOTTED 
NDPTS  -  ARRAY  OK  RCS  DATA  POINTS  IN  TIME  OR  RCSOTA 
NREC 1  -  POSITION  OF  FIRST  RECORD  READ  FROM  TAPE  -  TO  PRUVIDE 
T DENT IF  I C ATT  ON  INFORMATION  ON  THE  PLOTTED  OUTPUT 
HR  EC  2  -  POSITION  l)K  LAjT  REv-ORD  RcAD  KRfJM  TAPE  -  TO  PROVIDE 
IDENT  IFICATIUN  INFORMATION  AS  AbUVtl 
NPULSES  -  SMOOTHING  INIlKVAL  SELLCTLD  MY  USER  -  FUR  PLOT 
IDENTIFICATION 

NOUJ  -  OBJECT  NUMBER  -  FOR  PcOT  IDENTIFICATION 
FLG  -  A  FLAG  THAT  CONTROLS  AUTOMATIC  SYMETRICAl  SCALING  OF 
RCS  DATA  1Y  AXIS)  ABOUT  01. E.  IF  Flg- 1  *  SCALING  IS 
AUTOMATIC.  IF  NOT  1.  USER  MOST  SUPPLY  BEGINNING  AND 
ENDING  POINTS  OF  THL  Y  AXIS 

outputs: 

LINE  PLOT  OF  RADAR  CROSS  Sr.CTtPN  VS  ASPECT  ANGLE  OR  PROJECTED 
ASPECT  ANGLE  WITH  APPROPRIATE  I OuNT I F Y I NG  INFORMATION  ON 
CALS PA NS  CAlCOMP  PLOTTER 

SUBROUTINE  AND  FUNCTION  SUBPROGRAMS  REQUIRED  _ 

CALSPAN  PLOT  SUBROUTINES 

PLOT  -  BASIC  PLOT  SUB  USED  BY  HIGHER  LEVEL  PLOT  ROUTINES 
AXIS  -  DRAWS  A  OR  Y  AXIS  WITH  TIC  MARKS  AND  LABt.LS 
SYMBOL  -  DR*WS  AL->HA  CHARACTERS  QF  SPECIFIC  SIZE  AND  LOCATION  . 
NUMBER  -  DRAWS  NUMERIC  INFORMATION  OF  SPECIFIED  SIZE  AT 
A  SPECIFIED  LOCATION 

LINE  -  DRAWS  A  STRAIGHT  LINE  BETWEEN  POINTS  DEFINED  BY 
SUCCESSIVE  SEGMENTS  OF  X  AND  Y  ARRAYS 
IBM  SUPPLIED  FORTRAN  FUNCTIONS 

INT  -  RETURNS  INTEGER  PORTION  OF  FlGATING  PLUINT  NUMBER 
AMAX  1  -  RETURNS  MAXIMUM  VALUL  OF  ARG  OK  FLOATING  PT  NUMBERS 
AM  IN  I  -  RETURNS  MINIMUM  VALUE  OF  ARG  OF  FLOATING  PT  NUMBERS 

AUTHOR/ PROGRAMMER 
P.M  .MCMAHON 
APR  1979 
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